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1 Introduction 



Flavor-changing neutral current B meson decays based on the process b — > s£ + £~ 
(where I — e or \i) , such as B — >■ K*£ + £~ , provide promising probes for new physics. 
The decays proceed through loop diagrams such as those shown in Figure [U 




Figure 1: The electroweak penguin (left) and VT-box (right) diagrams responsible for 
B ->• K ( *H + £- decays. 



The theoretical treatment of b — > s£ + £~ transitions in the Standard Model (SM) 
follows an effective field theory approach in which the Hamiltonian is a sum of terms 
consisting CKM factors and Wilson coefficients that multiply operators formed from 
the light quark and lepton fields. The Wilson coefficients, obtained by integrating 
out the heavy particles, characterize the short- distance physics in these decays. New 
physics would modify the Wilson coefficients by providing new particles inside the 
loops and may modify the Hamiltonian by adding additional scalar or pseudoscalar 
terms. To account for QCD effects that mix the operators, so-called effective Wilson 
coefficients are defined. Measurements of b — > s£ + £~ decays probe the effective Wilson 
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coefficients Cf, Cf, and Cj£0 However, branching fraction measurements alone are 
not sufficient to exploit this opportunity, since they are generally in good agreement 
with SM theory predictions, even when measured versus q 2 (= mj e ) as shown in 
Figure [2] 




q 2 (GeV 2 /c 4 ) 

Figure 2: Partial branching fractions versus q 2 for B — > K£ + £~ (top) and B — > 
K*£+£- (bottom) for BABARpQ, Belief, CDF[3], and LHCb|H [5]. CDF and LHCb 
results are for modes only. The region between the magenta curves shows the 

SM range. [B] Yellow shading shows the charmonium exclusion region used by BABAR. 

More incisive observables are needed. Of particular interest is the lepton forward- 
backward asymmetry, Aps, which in the SM exhibits a predictable q 2 dependence that 
includes a zero-crossing point near q 2 = 4GeV 2 /c 4 . Non-SM processes may change 
the magnitudes and relative signs of the relevant Wilson coefficients, and induce 
relative complex phases between them, and in general may lead to large deviations 
of Afb versus q 2 from its SM expectation. 

1 Cj S represents the electromagnetic penguin; its magnitude is also probed by b — > s-f. Cg S 
and Ciq represent vector and axial- vector components, respectively, of the Z-penguin and W^-box 
diagrams. 
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In this note, BABAR's final (but still preliminary) angular analysis of the B — > 
K*£ + £~ decays is presented, based on a sample of 465 million BB pairs collected at 
the T(4S) at the SLAC PEP-II B-factory. This includes measurements of Fl, the 
fraction of longitudinal K* polarization, and lepton forward-backward asymmetry 
Afb in six q 2 bins (as defined in Tabled]). 



2 Data Analysis 

In this analysis B — > K£ + £~ and B — > K*£ + £~ events are reconstructed in nine dis- 
tinct submodes: B + ->• K + e + e~, B + K+fi+fj,-, B° K° s e + e~, B° ->• K^+p-; 
B + -± K+TT°e + e-, B+ K° s n + e + e-, B + ->• K° s tt + fM + ^ , B° K+ir-e+e', and 
5° — >■ K + n~ , with — >• 7T + 7r _ . Charge conjugation is implied here and 
throughout this note. The submode _B + — )■ K + 7T°{i + is excluded from this anal- 
ysis since Monte Carlo simulations indicated that it did not improve the results. 
Events where the dilepton pair originated from J /if) or ip(2S) decays are explicitly 
removed using selection criteria applied to m«. The major background sources are 
semi-leptonic B and D decays, which are suppressed by applying optimized cuts on 
bagged decision tree outputs. The bagged decision trees utilize event shape variables, 
vertex information, missing energy, and similar inputs to discriminate between sig- 
nal and background events. Training is carried out on Monte Carlo samples. An 
additional background from B + — > Dtt + , followed by D — >■ K*tc along with n /j, 
misidentification, is vetoed explicitly by rejecting K*n combinations near the D mass. 
After all selection criteria have been applied, the signal efficiency is typically about 
15%, although it varies by submode and g 2 -bin. 

After event selection, a sequence of three maximum likelihood fits is performed 
to determine signal yields in each g 2 -bin and to determine Fl and Afb- The first 
fit determines the number of signal events in each g 2 -bin. It is performed in three 



dimensions: tties, M(Kir), and C, where tties — y-^beam ~V*bi M(Ktt) is the mass 
of the K* candidate, and £ is a likelihood ratio formed from the output of the bagged 
decision trees used to separate signal events from other B meson decays. E^ eam and 
p* B are the beam energy and momentum of the B in the T(4S) frame (CM). 

The subsequent fits are performed for the signal enriched region defined by rriEs > 
5.27 GeV. In the second fit, a fourth dimension is added, cos 9k, where 6k is the angle 
between the K and the B in the K* rest frame. The yields determined from the first 
fit are fixed in the second fit, and the parameter Fl is determined based on the 
equation: 

It^tT = 1 Fl cos2 6k + I^ 1 - Fl ^ 1 - cos2 W 

1 a cos 6k 2 4 
In the final fit, Fl is fixed at the result of the previous fit, and a fifth dimension is 
added, cosd^, where 6g is the angle between the £ + (t~) and the B (B) in the dilepton 
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rest frame. The parameter Afb is determined from the fit based on the equation: 
1 c/r 3 3 

¥1 F = -F L (l-cos 2 ^) + -(l-F i )(l + cos 2 ^) + A FB cos^ (2) 

1 a cos 6 1 4 8 

The fit methodology was validated using both Monte Carlo test experiments and by 
fits to B — > K*J/ip and B — > K*ip(2S) modes in data for which prior measurements 
exist. 



3 Results 

Preliminary BABAR results for Fl and Afb versus q 2 are given in Tabled] While 
the dominant errors are statistical, systematic error estimates are included in Table [TJ 
Figure [3] shows these results along with results reported thus far by other experiments. 



q 2 bin (GeV 2 /c 4 



0.1 - 2.0 
2.0 - 4.3 
4.3 - 8.1 

10.1 - 12.9 

14.2 - 16.00 
> 16.00 



1.00 - 6.00 



F L 
=rrnr 



A 



FB 



0.23I£$ ± 0.04 

0.15lg;£ ± 0.04 
0.32 ±0.12 ±0.06 
0.40 ±0.12 ±0.06 
0.43±g;i§ ± 0.09 
0-55±8;if ± 0.03 



0.25l£$j ± 0.03 



0.14l^ ± 0.20 
0-40±gil ± 0.07 
0.15 ±0.16 ±0.08 
0.36l|J;$ ±0.10 
0.34±g;?| ± 0.07 
0.34l{$? ± 0.07 



0.17l^u ±0.07 



Table 1: Preliminary BABAR results for Fl and Afb versus q 2 . In addition Fl and 
A F b measurements are provided in a bin (1.00 - 6.00) that is sometimes used by 
theorists. The first error is statistical and the second is systematic. 



4 Conclusions 

Significant progress has been made in the study of b — > s£ + £~ related B decays. Thus 
far no significant deviations from the SM have been observed. Figure E] shows that 
the current results from all experiments are generally consistent with each other and 
with the SM expectations for Fl and Afb- However, for meaningful comparisons of 
these angular measurements with theory, much larger data sets are needed. LHCb 
measurements will improve as additional data is collected. Future super-B factory 
results will also be important since a larger number of final states will be accessible. 

This work was supported in part by Department of Energy contract DE-AC02- 
76SF00515. 
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Figure 3: Results for F L (top) and A FB (bottom) versus q 2 . BABAR preliminary 
results are indicated by the solid red squares. Results from other experiments are: 
Belle [2] blue stars, CDF [7] black triangles, and LHCb[5] solid black dots. Yellow 
shading indicates the charmonium exclusion region. The binned SM theory prediction 
[H] is shown in magenta. 
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